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Recently, it has been shown that the primordial black holes (PBHs) produced by near critical 
collapse in the expanding universe have a scaling mass relation similar to that of black holes produced 
in asymptotically flat spacetime. Distinct from PBHs formed with mass about the horizon mass 
(Type I), the PBHs with the scaling relation (Type II) can be created with a range of masses at a 
given formation time. In general, only the case in which the PBH formation is concentrated at one 
epoch has been considered. However, it is expected that PBH formation is possible over a broad 
range of epochs if the density fluctuation has a rather large amplitude and smooth scale dependence. 
In this paper, we study the PBH formation for both types assuming the power spectrum of double 
inflationary models in which the small scale fluctuations could have large amplitudes independent 
of the CMBR anisotropy. The mass spectrum of Type II PBHs is newly constructed without 
limiting the PBH formation period. The double inflationary power spectrum is assumed to be of 
double simple power-law which are smoothly connected. Under the assumed power spectrum, the 
accumulation of small PBHs formed at later times is important and the mass range is significantly 
broadened for both Types. The PBH mass spectra are far smoother than the observed MACHO 
spectrum due to our assumption of a smooth spectrum. In order to fit the observation, a more spiky 
spectrum is required. 
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I. INTRODUCTION 

The density fluctuations give an interesting way to 
form primordial black holes (PBHs) in the early universe. 
If the fluctuation amplitude of an overdense region is of 
order unity when the fluctuation enters into the cosmo- 
logical horizon, the overdense region can evolve into a 
black hole. Early studies found that the PBHs so formed 
have mass about the horizon mass at the formation time 
(Type I) |,|. The PBHs can cover the complete mass 
range of black holes considered nowadays: from a mass of 
about the Planck mass ~ 2x 10~ 5 g to the mass of a black 
hole in the galactic bulge ~ 10 s M . The Hawking evapo- 
ration g of PBHs with mass < 5 x 10 14 g can affect many 
early universe phenomena. In particular, PBHs with ini- 
tial mass ~ 5 x 10 14 g would presently be at the final stage 
of their evaporation and could be observed through their 
energetic particle emission. Although PBHs have not yet 
been detected, the upper limits on the PBH number den- 
sity have been extensively studied j|Jj|- When the fluc- 
tuations are normalized to the CMBR anisotropy, PBH 
formation strongly constrains the spectral index of the 
density fluctuations, n < 1.25 ^|J7). PBHs surviving to- 
day and the massive relics which may remain at the end 
of a PBH lifetime ||] are natural candidates for dark mat- 
ter. PBHs with mass ~ O.5M are also candidates for 
the observed massive compact halo objects (MACHOs) 
H|oJ. It has been proposed that PBHs with the MA- 
CHO mass can be produced in certain inflationary mod- 
els |pr|-|l^| or at the vanishing of the sound velocity at 
the cosmological quark-hadron phase transition |l5[ ]. If 
the MACHO PBHs form coalescing binaries, they could 
be detected by gravitational wave interferometers Jl6| . 

On the other hand, critical phenomena have been ex- 
hibited in the gravitational collapse which produces a 



black hole. This was originally found by Choptuik in the 
numerical study of the gravitational collapse of massless 
scalar fields [[17| . Especially, he found that for one param- 
eter family of field configurations with controlling param- 
eter p, the black hole mass scales as (p — p c ) lB near and 
above the critical parameter p c with a universal exponent 
7 S ~ 0.37. It has been shown subsequently that similar 
mass relations hold for other situations, such as gravita- 
tional waves [Q and radiation fluid [ pj| . Interestingly, 
all the studies have found that y s 0.37. The analytic 
explanation has been given by perturbative analysis [ p0| . 

Cosmological application of the critical phenomena 
was recently noticed by Niemeyer and Jedamzik in the 
study of the gravitational collapse of an overdense region 
in the expanding universe [ plp^ ]. They found that if 
the horizon crossing amplitude 5h of the overdense re- 
gion is near and above the critical amplitude Sjj w 0.7, 
then the black hole mass scales similarly, with j s « 0.36 
(Type II PBHs) [^2). The scaling relation can induce 
great changes in the mass spectrum because PBHs can 
form with a range of masses at a given time, while Type 
I PBHs have an initial mass of about the horizon mass at 
a given time. The mass spectrum of Type II PBHs has 
been studied for the case in which the PBHs form over 
a very short period [^l|,^3 24 1, such as occurs for a blue- 



shifted density perturbation spectrum or the spiky spec- 
trum associated with MACHO PBHs. Based on these 
mass spectra, the upper limits on the PBH number den- 
sity and the spectral index have been revised |24|]25j ] . 

With normalization of the density fluctuation spec- 
trum to the COBE CMBR anisotropy measurements, the 
PBHs can only form significantly under a simple power- 
law spectrum of the density fluctuations if the density 
fluctuation is largely blue-shifted with a fine-tuned spec- 
tral index |,g,||,f27|. In this case, the PBH formation is 
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concentrated at the initial time U when the fluctuation 
develops. Similarly, a narrow spiky fluctuation spectrum 
peaked at the MACHO scale is required for some MA- 
CHO PBH models [0,|l|,|l4]|. Hence, it has been usually 
assumed that PBHs are formed at only one epoch or with 
only one fluctuation scale. For Type 1 PBHs this assump- 
tion leads to the S- function type mass spectrum Q . For 
Type II PBHs, even if they are formed at one epoch, the 
PBHs can have a range of masses because of the scaling 
mass relation [^l| , p3tp4[ . 

With rather large fluctuation amplitudes and a smooth 
power spectrum, the PBH formation is possible over a 
range of epochs (or fluctuation scales) and the accumula- 
tion of PBHs formed at different times should be consid- 
ered. However, even in the models with large amplitude 
and smooth scale dependence j^], the accumulation ef- 
fect has previously been neglected and the mass spectrum 
obtained by simply taking a short duration of PBH for- 
mation. With an assumption of large amplitude on small 
scales, we will investigate the accumulation effect in the 
mass spectrum of PBHs obtained based on the Press- 
Schechter method p8| . For this, we newly construct the 
mass spectrum which is more adequate to describe Type 
II PBHs formed over a range of epochs. Related work on 
the mass spectrum was done in Ref. pj| using the excur- 
sion set formalism, but that reference investigated the 
validity of the mass spectrum obtained by limiting the 
PBH formation epoch. The large amplitude power spec- 
trum is natural in double inflationary models. In the 
double inflationary models, the small scale fluctuations 
could have large amplitude not constrained by the CMBR 
anisotropy, leading to enhanced formation of small PBHs. 
Though there are various double inflationary models, we 
will simply assume a double power-law spectrum. This 
kind of power spectrum has been studied in the dou- 
ble inflationary model of supergravity theory proposed 
to explain the MACHOs as PBHs Il2l. From the mass 
spectrum of PBHs of both types, it will be shown that 
the PBH mass spectrum has difficulty explaining the ob- 
served MACHO spectrum with the assumed power spec- 
trum of the density fluctuation. 

The rest of the paper is organized as follows. In Sec. 
II, the PBHs of both types are reviewed. The mass spec- 
trums are formulated in Sec. III. PBH formation under 
the double-inflationary power spectrum is given in Sec. 
IV. The paper closes with some concluding remarks in 
Sec. V. 



II. PBH FORMATION AND ITS MASS 

We only consider the PBH formation in a universe with 
a hard equation of state for which the sound velocity 
v s = ^7 with < 7 < I. Although we are interested 
in the radiation-dominated era (when 7 = 1/3), we re- 
gard 7 as a parameter for convenience. For an overdense 
region in the fluctuated universe to collapse, the fluctu- 



ation amplitude of the region should be large enough to 
overcome the Jeans pressure of the region. Conversely, 
it should not be so large that the region is decoupled 
from the background universe. These conditions can be 
written for the horizon crossing amplitude Sh as 

6 C H <S H <8% . (1) 

In addition, Sh should be larger than the critical value 
needed for the trapped surface formation. 

For Type I PBHs, it has been found from the study of 
the evolution of a spherical overdense region in the ex- 
panding universe that the bounds for the PBH formation, 
5 H and 5 H , are scale independent and should be of the 
order of 7 Q . It has been usually taken that S H =j and 
S H = 1. Also, it has been assumed that the trapped sur- 
face forms with these bounds at nearly the horizon time. 
Then the PBH mass can be approximated, neglecting the 
slight dependence on the fluctuation amplitude, as |2j[^] 

M BH (t) = 7 3 / 2 M H = j 3/2 M m (J-^J 

= 7 37/(l+37) M {l+l)/(l+3~)) M ^( 1+3 ^ (2) 

where Mh is the horizon mass, the subscript V repre- 
sents the quantity at ti, the time when the density fluc- 
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tuation develops, and M cx AIjj is the mass contained 
in the overdense region with comoving wavenumber k 
at ti. Note that the mass of PBHs formed at a given 
time is fixed as given in Eq. (2). So, all the quantities, 
Mbh, Mh, M , and k have one-to-one correspondences. 
Larger PBHs can form only if the larger fluctuations cross 
the horizon at later times. 

Distinct from Type I PBHs, Niemeyer and Jedamzik 
investigated the PBH formation by numerically simulat- 
ing the evolution of an overdense region in the expanding 
universe |pl| , p2| . Applying the three different configura- 
tions of the overdense region, they found that the critical 
phenomena also exist in the PBH formation and the re- 
sulting PBH mass scales as 

M BH = kM H (S H -S H )^ (3) 

where k ~ 3, 5 C H ~ 0.7, and the critical exponent, 
7 S ~ 0.36, is similar to that obtained in the asymptot- 
ically flat cases |22]] . In this case, the PBH formation 
mass is not unique at a given time any more but covers a 
range of masses, in principle, from zero to infinity. How- 
ever, we introduce the Planck mass Mpi as a minimum 
PBH mass. Though 5 H has not yet been determined for 
Type II PBHs, we take S H = 1 as the upper limit for 
PBH formation. Actually, S H is not crucial in this work. 
Then the possible mass range for a given time t is from 
Mpi to kM H {t){5 H -5 H P" ~ 2M H (t). Because the criti- 
cal amplitude 8 C H = 0.7 is larger than that taken for Type 
I PBHs 8 C H = 1/3, one may suspect that Type I PBHs 
always form more abundantly than Type II PBHs. How- 
ever this is not so because the mass shift to smaller PBHs 
should be considered for Type II PBHs. The details are 
given in Sec. IV. 
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III. MASS SPECTRUM OF PBHS 

The mass spectrum of PBHs can be constructed based 
on the Press-Shechter method To do this we start 
from the description of an overdense region at U . We only 
consider Gaussian fluctuations^] For an overdense region 
with mass scale M and size R, the smoothed density field 
5m is defined by 



J«W= / d 3 y6(x + y)W M (y) 



(4) 



where <5(x) = (p(x) — p)/p, p is the background energy 
density of the universe, and VTa/(x) is the smoothing 
window function of scale M. The dispersion crjy/ , the 
standard deviation of the density contrast of the regions 
with M, is given by 

°M = ^r<<&(*)) = TFT / l^| 2 ^ k 2 (M) (5) 



where Vw ~ R 3 denotes the effective volume filtered by 
Wmi and 5^ and Hk are the Fourier transforms of <$(x) 
and W^x), respectively. 

For Gaussian fluctuations, the probability that the 
region of size R has density contrast in the range of 
(Si + d5i,5i) is 



P(M, S,)dS 



1 



2-kom 



exp 



2^m 



(6) 



The fluctuation amplitudes at ti are related to the am- 
plitudes at the horizon crossing time as follows, 
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(7) 



Since the smoothed density field under the window 
with size R sees only the structures larger than i?, 



(8) 



F(M,5 C ) = / P(M,5i)dSi 



can be interpreted as the fraction of the overdense regions 
larger than R and with Si > 5 C . And the fraction of the 
regions with Si > S c in the range of (M, M + dM) is given 
by 



f(M, 8 c )dM 



OF 



dM 



(9) 



For Type I PBHs, from Eq. (1) and (7), the condition 
for PBH formation can be written for Si as [0 



*The PBH formation from non-Gaussian fluctuations have 
been studied in ||,|^,^3| 
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(10) 



From the PBH formation condition, the fraction of over- 
dense regions which will evolve to PBHs is fsndM = 
f(M,a)dM - f(M,0)dM ~ f(M,a)dM. Multiplying 
by Pi/M, this fraction can be converted to the number 
density. After a change of variables, the initial mass spec- 
trum of PBHs in the range of (Mbh, Mbh + dMsH) is 
given by |fl 



n BH (M B H)dM BH 
'M B h do H 
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BH 



(11) 



Here a factor of 2 is included in accordance with the 
Press-Schechter prescription for normalization p8[ . The 
number density of all the PBHs at time t is then given 
by 



0BH(t) = 



a(t) 
a (ft) 



Mb 



Mb hi 



n B H(M BH )dM B H (12) 



where a(t) is the cosmic scale factor, Mbhi — 
j 3 / 2 M H (t), and M BHX = M m if M.(t) ~ 
10 9 (t/lsec) 1/3 g < M m where M»(i) is the PBH mass 
whose lifetime is t, else Mbhi — M*(t). If PBHs leave 
stable massive relics at the end of their evaporation life- 
times, then Mbhi = M Hi- 

If the power spectrum can be approximately treated as 
a ^-function type as in the blue-shift perturbation with 
the CMBR normalization or the spiky spectrum in some 
inflationary models, it is usual to assume that the PBH 
formation occurs at only one epoch. In this case, the 
PBHs form with only one mass for Type I PBHs, and 
the initial number density is approximately expressed by 



Pt(M B H) = oh exp 



2a H 



(13) 



For Type II PBHs, even if the power spectrum resembles 
a (5-function, the PBHs can form with a range of masses 
from about Mpi to 2Mn(t). Their mass spectrum can be 
obtained from the probability distribution of Si (or Sh), 
Eq. (6), as follows, 



tibh (Mbh ) dMsH = ■ 
(M BH \ lhs 

x exp 



i 



Pi ( M H 



2tt<thJs Mhi \Mni 



-3/2 



S 



H 



dM 



BH 



2a H 2 Mbh 



where, from Eq.(3), 



oh = o H - 



(M BH V /7s 
\kM H ) 



(14) 



(15) 
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The excursion set formalism of [£3| validated Eq. (14) for 
the blue-shifted fluctuation spectrum with small CMBR 
anisotropy and the spiky spectrum for MACHO PBHs 
p3| . However, the above expressions Eq. (13) and (14) 
are inappropriate to describe the accumulation of PBHs 
formed at later times. In particular, for Type II PBHs, 
the accumulation of newly formed PBHs contributes to 
all masses below ~ 2Mjj(t), so that the shape of the mass 
spectrum is time dependent while that of Type I PBHs, 
Eq. (11), does not depend on time. For Type I PBHs, 
the newly formed PBHs just add to the total number 
density at masses higher than the older PBHs. 

One of the main purpose of this work is to formulate 
the mass spectrum for Type II PBHs which describes 
the accumulation effect more properly. For this we start 
from the interpretation of the probability distribution, 
Eq.(6). As mentioned before, the smoothed field can see 
the structures larger than the window size. So, Eq.(6) 
can be interpreted as the fraction of the overdense regions 
which are larger than the scale M and have amplitude in 
the range of (Si, Si + dSi). The contribution from larger 
scales can be eliminated by differentiation. Then the 
fraction of the regions between (M, M + dM) and (Si, Si + 
dSi) is given by 

f(M, Si)d5idM = - dP ^' 5l) dSidM . (16) 

Applying a change of variables (M — > Mh, Si — > Mbh) 
and the Press-Schechter prescription, we can find the 
mass spectrum for the PBHs. The number density of 
PBHs in the ranges of (Mbh, Mbh + dMsn) and 
(Mh, M h + dM H ) can be derived to be 



'2 1 

h B HdM H dM B H = -\l — 



M 



H 



-3/2 



irjs M m \M m 
' ' < \ M//, da H 1 



/ Mbh\ ( Sh_ _ 
\kM H ) \<jh 2 



ojj dM H o h 



x exp 



S H 2 \ dM H dM BH 
2a% 



M H Mbh 

The mass spectrum at time t is then given by 



(17) 



n B H(M B H,t)dM B H 



M h 



n BH dM H I dMBH (18) 



where M m = M BH /[k(S d H - S L H )^} and M H2 = M H (t). 
The number density of PBHs at t is 



0BH(t) = 



a(U) 



3 pM BH2 
J M B hi 



n BH (M b n)dM b h (19) 



where M BH 2 = kM H (t)(5 d H - 5%)^ and M B hi = M Pl if 
M*(i) < Mp h otherwise M B hi = M*(t). If PBHs leave 
stable relics, then we can take Mb hi = M re n c . 



IV. PBHS UNDER THE DOUBLE 
INFLATIONARY POWER SPECTRUM 

With the COBE normalization on the CMBR 
anisotropy and a simple power-law spectrum, the PBHs 
can form significantly only if the spectrum is strongly 
blue-shifted. The PBH formation is then usually as- 
sumed to occur only at one epoch and the accumulation 
of black holes formed at later times is neglected. How- 
ever if the power spectrum can have larger amplitude, 
it is expected that the PBH formation is possible over a 
rather broad period of times. For practical purposes, one 
would require the large amplitude to occur on scales suffi- 
ciently small as to be irrelevant to the CMBR anisotropy 
while the fluctuation amplitudes on large scales are con- 
strained by the CMBR anisotropy 5h ~ 10~ 5 . This kind 
of power spectrum can be generated by double inflation- 
ary models |l^[l4|,|3l]|32| . In particular, the double infla- 
tion model with one scalar field has been studied in Rcfs. 
Jl3| , pTt . In the double inflationary models, the first in- 
flation gives the density fluctuation which is responsible 
for the CMBR anisotropy on large scales. The second 
inflation generates the smaller scale density fluctuations 
which are irrelevant to the CMBR anisotropy. If the 
small scale fluctuations have a rather large amplitude, 
then the PBH formation can be enhanced. Phenomeno- 
logically, the large amplitude on small scales in the double 
inflation model could resolve the problems of the stan- 
dard cold dark matter scenario plj . It has been also 
claimed that MACHO PBHs could be produced in mod- 
els with large amplitude at the MACHO scale ~ lM r . 
such as the double inflation from supergravity theory 
and chaotic new inflationary potential Jl^,|lJ] . The dou- 
ble inflation models with two scalar fields might suffer 
from the initial condition problem if our universe origi- 
nated from an inflationary domain of Planckian density 
p3[ . However, the study on this problem is not the sub- 
ject of this work and we assume that the double inflation 
is realized in a proper way. 

We consider only the fluctuation scales larger than 
Mm, the horizon scale at the reheating time of the sec- 
ond inflation (= ti). The PBHs from the fluctuations 
smaller than Mm, if an Y, will be diluted away by the 
second inflation. Also, for simplicity, we assume that the 
power spectrum larger than Mm is a combination of two 
simple power-law spectra which are given by 



(Th(Mh) = (J Hi 



M. 



H 



Mm 



(l-n)/4 



for M H < M C H 



f M \ ( 1 ~™° bs )/ 4 
<W tt 5 - 1 fOT M H -»M C H (20) 



\M_ 



HO 



where the present values Mho ~ 10 5b g, ano — 10 
§ and n obs = 1.2 ± 0.3 from the COBE data @, 
and the breaking scale M^ is well below 10 15 M Q which 
is the horizon mass at the matter-radiation equal time 
plf. The two simple power-law spectra are assumed to 
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be connected in such a way that an decreases mono- 
tonically to the second spectrum iorM B > Mfj. Also 
the change of an is assumed to be smooth satisfying 
(Mh /cni)(d<JH /dMff) <C 1. If the spectra were con- 
nected sharply or discontinuously at Mjj , PBHs or other 
structures with the mass of about Mfj might be pro- 
duced significantly. If Mjj ~ 10 33 g, such PBHs could be 
related to the MACHOs. However, with the above as- 
sumptions, we will not consider such PBH formation in 
this work. Also, by assuming n b s — 1, the PBHs from 
the fluctuation larger than Mfj are neglected. 

Then the mass spectrum of Type I PBHs with the 
assumed power spectrum is given by 

our 71 + 3 ft 1 Pi (MbhY^ 

n BH (M BH )dM BH = ^-J 

1 / 7 2 \dM BH 

x ^ ex V^J^T • (21) 

Fig. 1 shows the mass spectrums for various parameters. 
The PBH mass starts from M BHi = j 3 / 2 M H (U). For an 
n = 1 scale invariant spectrum, the mass spectrum has 
no cutoffs and scales as n BB oc M BH ~ 5/2 - For n < 1 
(n > 1), there are exponential cutoffs at small (large) 
scales. However, as the fluctuation amplitude grows or 
as n approaches 1, the cutoffs develop more slowly and 
the mass range is significantly broadened. 
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FIG. 1. The mass spectrum of Type I PBHs today for 
T RH = 10 10 GeV and M C H = 1O 5 M ; (a) for n > 1, I: n = 1.0, 
o Hi = 0.6, II: n = 1.0, o m = 0.06, III: n = 1.2, a Hi = 0.6, 
and IV: n = 1.2, om = 0.06, (b) for n < 1, I: n = 0.4, 
cth(Mmacho) = 0.6, II: n = 0.4, <th{M M acho) = 0.06, 
III: n = 0.8, <th{M M acho) = 0.6, and IV: n = 0.8, 
<jh{Mmacho) = 0.06. 



For Type II PBHs, the mass spectrum is more com- 
plicated. With the assumed double power spectrum, Eq. 
(17) becomes 
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FIG. 2. The mass spectrum of Type II PBHs today with 
the same parameter set of Fig. 1. 

Fig. 2 shows the mass spectrum ti B h{M B h 1 to)- The 
PBH mass spectrum extends down to M Pi . There is 
an extra condition for Type II PBHs: requiring posi- 
tive number density requires Sfj > a\. This condition 
causes a rapid drop in the mass spectrum at small scales 
for n > 1 or at large scales for n < 1. The drop can 
be seen in III of Fig. 2a. The mass range is also sig- 
nificantly broadened when the fluctuation amplitude is 

large or n approaches 1. Since 5 B is a function of M BB , 

s 2 

the role of the exponential term exp(— t^t) is very diffcr- 

2 

ent to that of cxp(— -^r) in the Type I mass spectrum 

which leads to cutoffs in the mass spectrum. Interest- 
ingly, the mass spectrum for Type II PBHs scales in the 

— 5/2 

same manner as for Type I PBHs, n BB oc M B ^ , in 
spite of the enhanced formation of small PBHs. This is 

s 2 

mainly because exp(— ^r) is a function of (M BB / M B ) 
for n = 1 and the main contribution to / h BB dM H comes 
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from Mbh ~ Mjj- For n > 1, the exponential cut- 
off develops similarly at large PBH mass and there is 
an Mbh 5 ^ 2 decrease at small PBH mass. However for 
ri < I there is no exponential cutoff at small mass be- 
cause of the accumulation of small PBHs formed at later 
times. Instead, the mass spectrum at small mass scales 

as ubh oc Mbh 1b ^^° ■ 

For both Types, the dependence on the reheating time 
only appears in the mass spectrum for n > 1. The mass 
spectrum shifts to small mass as the reheating time de- 
creases. Unfortunately, we now have an extra parameter 
M H . The breaking scale Mfj must be determined, in ad- 
dition to the shape of the power spectrum, to estimate 
the number density (3bh or the density fraction £Ibh 
of PBHs. This could be done for more concrete mod- 
els. Although (3bh and Qbh are not determined here, 
it is interesting to compare the mass spectrum of PBHs 
with that of MACHOs. From the data fit of 8 MACHO 
events, the fraction of MACHOs between Mmacho and 
Mmacho + (IMmacho is found to be 

/m AC HO dMM AC HO oc M M 3 ^ CHO dMMACHO (23) 

for Mmacho > 0.3M Q pc| |. This fraction corresponds 
to a PBH mass spectrum which is far steeper than the 
spectra for n < 1 shown in Fig. Ib and Fig. 2b. In Ref. 
p2| , it was claimed that the second inflation with an n < 
I spectrum could produce PBHs peaked at the MACHO 
mass. However, they did not consider the accumulation 
effect. It seems difficult to generate the current MACHO 
spectrum with the assumed two simple power-law spectra 
unless the spectra are connected in a proper way near 
Mmacho- The MACHO spectrum can arise for some 
n > 1 spectra, but the universe would be overdosed by 
the PBHs. 



V. CONCLUDING REMARKS 

In this paper, we study the formation of PBHs of two 
different types, the PBHs with horizon mass (Type I) and 
the PBHs with scaling mass relation (Type II) , under the 
double inflationary power spectrum. The assumed power 
spectrum is a double simple power-law spectrum such 
that the small scale power spectrum gives PBH forma- 
tion and large scale power spectrum generates the CMBR 
anisotropy. With a rather large fluctuation amplitude, 
the PBH formation occurs over a range of epochs and it 
is important to consider the accumulation of smaller mass 
PBHs formed at later times. For both types, the PBH 
mass ranges are significantly broadened. As expected, 
there are cutoffs which develop slowly in the mass spec- 
trum for Type I PBHs. From the newly constructed mass 
spectrum for Type II PBHs, it is found that the mass 
spectrum scales in the same manner as Type I PBHs for 
n = 1. Due to the enhanced formation of small PBHs 
for Type II PBHs, the mass spectrum for n < 1 does not 
develop an exponential cutoff at small mass but scale as 
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bh ^' 7s . This may give an unique feature of PBHs 
in the history of the universe. 

With the assumed power spectrum, the PBH mass 
spectrum has difficulty explaining the current MACHO 
spectrum. In the present work, however, we did not con- 
sider the spiky spectrum [|ll"||ll^ , |l4f| which has more rapid 
scale dependence. Since the resulting mass spectrum will 
also have a rapid scale dependence, one could explain the 
MACHO spectrum by controlling the model parameters. 
However, the Press-Schechter method could not be appli- 
cable if the <7h grows rapidly as can be seen for a simple 
power-law with n < —3. So, our mass spectrum could not 
be applicable to the spiky spectrum. This problematic 
issue will be studied in detail later. 

The gamma-rays from evaporating PBHs can gener- 
ate a diffuse 7-ray background (DGB) in the universe 
||,[35-37|. However, it is difficult to generate the complete 
observed DGB in the 7-ray energy range 0.8MeV < E < 
lOOGeV [B8,p9| with the PBHs from a simple power-law 
spectrum |7, [40]. For E > fOOMcV, the PBH 7-ray spec- 
trum falls off E~ 3 , faster than the observed DGB spec- 
trum, independent of the mass spectrum of PBHs |?5 36 
(or E~ 4 if the PBH emission can form a photosphere 
around the PBH pi]]). On the other hand, the 7-ray spec- 
trum for E < lOOMeV depends on the mass spectrum of 
PBHs in the range of 2 x 10 13 g < M BH < 5 x 10 14 g @. 
The double inflation model could be tuned to form suf- 
ficient PBHs to explain the observed E < lOOMeV DGB 
spectrum. 
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